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abstract 

We present a description of the data reduction and mLpmaking pipeline used for the 2008 observing season 
of the Atacama Cosmology Telescope (AC l). The d*u prevented here at MS Clti represent 12% of the 90TB 
collected by AC!' from 2007 to 2010. In 2008 we observed for 136 days* producing a total of 1423 h of dala 
( 1 1 TB for tiic 148 till/ band only), wilh a daily average or 10.5 h of observation. Prom iJhesc, 1085 h were 
devoted to a 8 50 deg* stripe 1 1 L2 h by 9? I j ix rile red on ll tketinahon of *52 P7, while 175 h were devoted to a 
280dqg ? stripe (4.5 h by 4TS) centered at the celestial equator. We discus spurred of utatistica] and systematic 
noise, calibration. telescope pointing, and datii scSeciioti. Out of 1260 survey hours and 1024 detectom per 
array, 816h and 593 effective detectors remain after dtirji selection for this frequency hand, yielding a ^8% 
survey efficiency. The total sensitivity in 2tK)3, determined from die noise level between 5 Hz and 20 Hz in die 
time-ordered dala stream (TOD). is32^Ky^ in CM@ units. Atmospheric brighincss llucuiaiions constiturc the 
main contaminant in the data and dominate the detector noise covariance at low frequencies in the TOD. The 
maps were made by solving the least squares problem using the Preconditioned Conjugate Gradient method, 
incorporating the details of the detector and noise correlations. Cras^correlatior: with WMAP nky [flapp., 
a* well as an is lysis front Simulations, reveal that our maps are unbiased at niuthpofe* t > 300. This paper 
accompanies the public release of the 148 GHz tombem stnpe maps from 200S. fhe techniques described here 
wft] be applied, lu future maps and dfua releases. 

Subject headings: Microwave Telescopes, CM13 Observations 
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I INTRODUCTION 

Over the past I wo decades precision measurements nf ihe 
Coy m; Microwave Background (CMB) have led to remark- 
able advances in nur understanding of cosmology, Combined 
with other observations, they have produced constraint* on 
models of the Universe to percent level accuracy (c.g.. Kjj- 
mal$u 21 al. (20! 1 )). Primary CMB anisotropy has been mea- 
sured to cosmic variance precision by WMAP up to multi-' 
poles of approximately 300 (Larson ct al 201 1 ). OnservaLions 
at finer angular scales (Friedman et a], 2000; Reichordt el al. 
200%, a; Vcncziani et al. 2D09; Kessler ct ul, 2009; Sicvtrs 
ct ill. 2009; Sharp ct aL 2010; ShirokofT et al. 2011; Das ct al. 
201 lb, Rcj chart* ct al 201 I f cmTcspi ru ding :n ihc dumping 
sea le of the anisotropiesi, have ted lq even tighter parmme 
ter constraints, which are improving- with each new data set 
Additionally. CMB meisurerneras at smaller angular scales 
-£C set \ihsc ti> :tu* ccnmbution oi point winch help 

to reveal the nature of early galaxie* or active galactic nu- 
clei (Vieira et al. 2010; Marriage ct al. 201 la), in (tie the: 
mal Sunyjuv-Zel'dDvich (SZ) effect, which can provide ;ndr 
render r CdnsirdinUs on cosmological parameters (Sun vac v & 
Zddorich 1970; Marriage ct al, 2011b: Schgal ct al. 2011, 
Planck Collaboration VtH 201 i; RciehardE ct al 20 1 2k atKl 
to graviiatiunul tensing effects on the CMB (Sctj uk 19%; Zul 
daniaga & Seljak 1999; Da&etal, 201 la; Shcrwin et al. 201 1), 
Fhc Atacama Cosmology Telescope (ACT) is located tit 
22 & 57'35"$, 67^47' 1 3" W on Cwro Tbco at an altitude of 
.“1200 m in the Atacama Desert of northern Chile. Its main pur- 
poses ..re to map the millimeter wave sky at arc minute scutes, 
sampling muhipolcs up to / ^ I0 4 and detecting and churac 
leri/ing foregrounds, including gfllaxy clusters through their 
SZ signature, and millimeter galaxies. 

Between 2007 and 20 U, ACT was equipped with the 
Millimeter Bolqmetric Array Camera (MBAC). which ob- 
served simultaneously in three bands: 148 GHz, 218(111/ .-md 
277 GHz. The bonds were chosen to avoid major atmospheric 
emission lines and to sample the SZ decrement, null and in- 
crement. Each band had a dedicated set of optic* and a detec- 
tor array composed of 1024 pop-up TES bolometers (Etenford 
el al, 2009) coupled to the optical signal, called "live deloe 
tors," p|u* J2 “dark detectors,” whreh were not coupled to the 
sky. Per details about the instrument see Swetz et al. (2011 1 
Since first ligtu on October 23, 2007, the telescope has 
had four observing sca.son.s producing more than 90 IB of 
data Here we describe the II TE3 of data taken at 1443 GHz in 
2008. including ihe techniques used for data characterization, 
ns well ns the data reduction process ured to obtain Eire filial 
maps, ‘Hk 218 GHz ami 277 GHz data wi!t tie described in a 
tetei' paper. However, a similar data reduction pipeline is used 
to analyze the full dataset. 

The maps obtained from these data hove an angular resrHu- 
lion of 1 .37' (Hindus cl al. 2010) and a noise level that ranges 
between 25 and 50jrK aremin. The calibration of the data 
to WMAP it discussed in Hajion et al. (201 1). The power 
spectrum of the maps is presented in Fowler et al, (20!Ch and 
Dh et al. •: 201 lb), with corresponding constraims cm con 
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TABLE i 

SCANPAftAMGTttlf 


Season 

2007 


JliOV-JifiL: 

Elevation 

3Q*S 

50*5 

50° 5 

jKian Width 

*i?b 

TP 1 !) 


Etnod 

1-9 4* 

1(I.?H 

lUJs 

Speed 

IJJtktt/l 

1 5des/a 

1 .5 tiep/a. 

Max. Accel, 

B.ttkg/j 1 

3 lileu/v 1 

3, J cIiuk/’i 2 

Dils hie Ira |?ih 

15 min 

35 flrJn 

10 till ii 


m otas i cal parameters in Dunk Icy cl ul. (2011). Cluster de- 
tections through their SZ signature arc presented by Hindu 
et al. (20| 0i and Marriage et al. (201 lb), white erftrapdaclic 
source detections are given in Marriage ct al (20 II a), Multi- 
wavekngth follow up for these clusters k presented in Menan- 
tcau et a). (2010) os well os the discovery a massive cluster at 
z - 0.87, (Mmonieau et al. 2012^ EL Gonfcjk Their cwmo- 
logiral interpretation is discussed in Schgal ct al. (201 1) The 
first di retd detection of gravitational tending oi tine microwave 
background was made using these mops in Dak et al. (201 la). 
This in turn demonstrates for the first time liuu microwave 
background data on their own favor cosmo logie-s with an ac- 
celerating expansion (Sherwincl al. 201 I h 

This paper is organized as follows, A summary of ob- 
servations is given in Section 2, In Section 3 we provide 
lmck ground for understanding the sky signal hs recorded in 
ihc time-ordered data stream (TOD), In Sections 4, 5 and 6 
respectively, we characterize the atmospheric, detector and 
systematic noise found in the TDD Datu calibration inlo 
imiLs of sty temperarm'e fluctuations. (STcmti, is described m 
Section 7. Section 8 do scribes the pointing solution, while 
Secdon 9 explains the detector time cunHtant dcicr ruination 
rneihod. fJatd ic lection is described in Section |0 pro rid 
■ng final statistics on the amount of data used for making the 
maps. The map making method is discussed in Section 1 1 . We 
conclude in Section 12 and present a stqv-by step summary of 
the data pipeline from raw data U> nuipv. Appendix A provides 
fiuirther details about finding and removing condoled modes 
from ihc data. 

This paper accompanies the public release of the data 
through the NASA's LAMBDA site. 12 

2. OBSERVATIONS 

AGT ob.xcrvcfi the ?ky by scanning ihc telescope in azimuth 
ut li constant elevation of 50 n 5 as the sky moves across the 
Ml Id of view in time, resulting in a stripe shaped obtteTVO* 
lion area. With this scan ^iraEegy, tItc instrument observes 
through li constani air mass, the cryogenic, remnin stable, the 
telescope's shape remains, constant, and ihe local environment 
ind insrmmental offsets are sampled in ,i comti.sLcnt way riie 
time constants of the detectors, together with mechanical fuc 
tors, limit the scan speed and turnaround acceleration. Table 1 
gives a summary of the scan parameters used in all &casuns. 
The lower acceleration introduced m 20(18 was needed to re- 
duce vibrations in the optics (see mote details in frfi). 

The obitervatirms are repealed at complementary cemral a^- 
itnudi angles to caprure both rising and setting skies;. This 
cross-linking technique helps minimize systc malic clfccis due 
to ihc --cun in the map making process and improves the deict- 
mmaiitm of CMB modes parallel to the scan direction. 

I! anp;^/l govi'prod'^c'u bcl 
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t if. l EUamptc TOD fnu» one iktcMU* «l 14KCi Hu from October 
2l„ lOft-i I'tpm wh ■ flood otHervinjg nighi with a PWV of 0.22 men. 
The tfu* drift it dnflunaJtidty changr* m the atmosphere brifhlticM. 
The biftifrcquency otHie ti dominated by deiettof unite. Dojis <ut 
mK in ,:MK equivalent ultiLi me I4f(£5lii: The plot display* 1,6 k 
ll/ unpin at nn inlcrvul of 2 5 y HO" J wond* The teieMropc 
m- mining whde these dm* were taken. 


Dach detector is campled it a rai£ of 398.72 Hz. nmJ Ihu 
did* me stored in 1 5- minute long deu files. They nrr then 
mni'gct. with the next of the housekeeping dull*. which in- 
clude ihe azimuth nrcd elevation encoder rc.iiJiii^ urn! lirnc 
of day. The dnui Lie comprised to one-third of their Ortgi- 
n.ii i: for 'iiiM age, munp t lossless eoniprrvninri iilgoMihrn 
{SLIMj. 11 A sunpte TOD b shown in Figure I far b night 
with goad weather conditions, meaning that pranptiahte wa 
Ujf vapor |PWV) remained bdow I mm The P W V i* a good 
indicator of 148 GK/ opacity and overall dam quality. 

Despite wime variation oscr the season due to changes In 
untrue und sumri times a typical day of ohservatiafUi at dir 
Kite wit ax follows. The cryogenic system wo* recycled ev- 
ery lUv before providing roughly J4 hottr> of cfcwTving tunc 
The cr/ogcnk cycle began around I TOO in the morning to 
nl time and lasted for 9.5 hours, so that MB AC was cold by 
about 20; 30, At 20:40, warm up movement* fur tin: motor* 
and gran were run. and at 20:50 the detector* were biased 
and liie lirst del ecu ir calibnition data were obtained (see |7 I } 
The observations Mimed at 2 1 :00. usually by scanning the lift 
trig sky. Around 2:30 mane detector calibration data were col- 
tooled and die *cun was shifted to Ihe west to observe dm ume 
legion that wu previously observed white rising Observing 
m die v c vt frA*. an jdditioncl .id vantage from a hardware safety 
perspective: in the event of a telescope failure, MBAC would 
he pouting *w»y from ihe Sun at its rising- Qbteivaljunsi tiof- 
maJly ended around H>4I). nearly 2 hours after tunnte abive 
the mountain*. ** At das point, final detector calibration data 
w ere liken, before the telescope was sent to it* hone pmu- 
Uun mJ the cycle was restarted When ihe otaervahtt region 
of ihe sky contained i planet, it was normally vanned every 
other night for calibration and beam measure menu, This was 

" 1 hut funder item ! Mr aIURlilT,*i - en. r :* 1 cs 1 3* nr 

114 ftof poinlinf iu4iii[7 muni we ended up wnfq dm ibnflfd 
k-u ihui oor Ihh; after itnrar 


TABLE 2 

CtaS£*v^TKIK tUKKABT FOt Et ASOKt 2000 {*1 14J GHCi 


Smhb Ite 'Mill Sin IU;Min, Mn> Mel ! Jiy 1 ) IT 


Wibrirw Sfiipr 




20M (-5T*!&-«?n 

I2CT<3*. 7^33*) 

050 

I0tl5 

t-quahirini Stripe 




21KH (-2=] 2, 2, S 34J 

(Mfai 1 *. 14"^) 

200 

ns 


“ T-Jlfll Imyni hefoii; data iCfetOn. 



Fm. 2.— Histogram showing the PWV during the 2008 sea- 
son. steu&jrcd a: the zenith nf the AfiKWU Pittilinder EiperLmeni 
(APFJCJ facility. The median value U 0.49 mm. 


not done every single night to avoid producing a gup in the 
CAfR map. AH of the operations listed above woe automated 
and could be pcifonncd irmotcly. 

The 2008 season began on AoguA 1 1 and ended on Decenv 
ber 34. with i itnal of 1 36 available my his and 124 nighis with 
■uccctafuJ observitoon^ tevultrag in 1423 hours and 25.7 TB 
of total dal* for all three freqLcncica, Prom the total observed 
nights, four had bad weather conditions, leaving 120 nights 
with usable science observations. The overall calendar lime 
efficiency, including day time, was 44% Pie median PWV 
during observations across the muon, measured at zenith of 
the Atacama Pathfinder Lx pen mem (AFFIX) facility (Glisten 
et al, 3006^, was (1.49 mm Given that the ACT site is abcHt 
1-Wlrn hug her (han the APFX vite. we use • cottefion factor 
PWVact - O.SRPWV^nrx n> c*.litrwtc the PWV at the ACT 
sale Figure 2 shows a hisingram ol the PWV during the sea 
son. 

Observations were made in tvro areas of the sky : The equa- 
torial snipe, centered at a declination at 0* (60* and 3(XF io 
uimtiift at 50 5 elevation i. and the southern stripe, centered 
n ■ declination of -53 d f 1 50“ and ’lO^iumuth at 50f5 d- 
evation). covering a wide range m right ascension Table 2 
lists the boundaries of the observation areas and the immBmj 
of hours available in each area before data selection. 

J r SIGMA! PHQPIiHTlES 

In tMi section wc characterize Ihe expected TDD signal 
from point sources and extended source*, for comparison to 
Ihe random und systematic noise dtNcnhed Idler in Sections 5 
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and 6. 


3.1. Pm nt Sources 


Fra "hji purpose*, wc will define a point source as any ob 
jeci with an angular size comparable to or smaller than the 
beam of The telescope. The beam full- width ai half max- 
imum (IAVHM), 8 : 3 , i\ L r 37 for 148 GHj.. Planets are the 
most important point sources as their high surface brightness 
makes them useful fur both calibration and beam measure 
menus. Distant galaxies can also be approximated as point 
sources smd are helplut for the pointing solution. 

As the telescope scans and the sky rotates, a single detector 
traces a zigzag on the sky* sampling the point source every 
time it intersects with the beam. The number of rimes a point 
source appears depends cm the scan period, the beam size and 
the central scan azimuth. A point source appears as a succes- 
sion of "blips" in the time stream. The shape of the*c "blips" 
is a slice through the telescope point spread function and de 
pends upon (lie angular separation between the center of the 
beam and (he location of the source The angular speed of 
the scan on the sky is given by 8 *= v^coslSOf 5), The 2008 
scan speed of v iafi - I ?5fs. implies an singular sky spece of 
57^2/s, Assuming a □uussiau-like beam of equivalent width 
and neglecting the transit speed of the source, the 3 dll cutoff 
frequency is 


/ .WE - 


2\n{2}0 

w8 l/z 


(I) 


This mean* that Tlvc contribution from point sources to ifae 
TOD is limited to frequencies below 1 8 Mv., as can be seen 
in Figure 3, 

C riven our scan speed, nominal elevation, sky rotation and 
sampling rate, wc expect nearly 10 samples per beam on the 
sky at 148 GHz. 


3,2. Exttnekd Soutt&s 

Hung CMB simulations we can csti mute the expected TOD 
response to extended sky features. Averaging the power spec- 
tra from many such synthetic TODs, we can estimate the 
CMB power spectrum m TOD spate . Figure 3 shows The av- 
erage: power qjcctrunL of u simulated UK GHz observation of 
the CMB sky and the average data power spectrum from one 
15-minute file. Because of Silk dumping, the CMB has a ml- 
atively sharp cutoff in ll$ characteristic size: unlike the point 
source observations shown in The figure, the CMB has Jilile 
si gnu! at frequencies higher than IflHz. 

Given ACTs angular scan .speed, the TOD frequency asso- 
ciated with a iiky feature of angular size 8j is 

£ 

(2> 

where the factor of 2 in the denominator comes from consid- 
ering that the angular wale is the size of a positive or nega- 
tive temperature bump, which is half or a wavelength on the 
sky. Given that multipole moment £ relates to angular wade 
us f ~ ir/fy for scale* x trial I compared to the full sky, m ap- 
proximate convention between multipole and TOD frequency 

U 

f ^frr/bp„ 2*fTQ o m 

0 v to cos(5Q?5J 

Note that this relation is only a rule of thumb: the actual 
mapping of TDD frequency into multipole space depends on 


the specific scan strategy This relation shows rhai the CMB 
power spectrum in the TOD can be shifted in frequency by 
changing the scan rale. As a reference* the TOD frequency 
corresponding to £=3000 was 7 9 Hz in the 2008 season 


i 

]fl> m 1 kt Lf^ 



FrEOLency 

Fra. 3 — Powcc spectra uf higna^ rn the seasen compared to 
the hi* la- The solid Curve is the average power spec Hum from 737 
live detector TODs from one I5~miaute 148 OH* data file, during 
which Lhc PWV was 0.5 mm. The rise at low frequencies js the noise 
contribution from die atmosphere. The dashed Curve stows the siro- 
ii.3j.Urd. C!MH tignsl Ihc oscilbEjotls in die simulated OMR sjgnpl 
are dm? to tmhanoed power at htrmonics of the scan freque-ticy. The 
dol-dashed curve estimates die point source contribution (consider- 
ing. only one hit at the beam center* assuming a Gaussian beam, for 
comparison, the thin and thick dotted curves show the expected re- 
sponse for 218GHz and 277GHz, with beam sizes of I fO0 and 0!9l 
respective fy, TTic umpiLtude of die point source power spectra com®- 
npundi to approsimatety the signal of Saturn. 


4. ATMOSPHERIC EMISSION 

The atmosphere is a strong emitter and absorber at Uic 
bands of i merest, chiefly due iu the excitation of the vibra- 
tion at and rotational mtidcti of water vapor. For this reason, 
rhe PWV is stmugiy correlated to the level of optical load 
ing on the detectors. The Atmospheric Transmission at Mi- 
crowaves (ATM) model (Pardo et al. 2001) provides an esti- 
mate of The loading as a function of the PWV level. At the 
median value of 0,40 mm during the 2008 observations, the 
loading is approximately 0.5 pW, wtiich corresponds to an 
equivalent Rayleigh- Jeans temperature of 6.4 K at the nom 
inul elevation. 

During the reason, the median atmospheric temperature 
drift uver 1 5-ininute observations was 0.22 K i Ray Leigh leans 
equivalent enits as measured hy MB AC). with lower and up- 
per qualities at 0. ID K and 0,43 K respectively. 

Turbo knee induces a spatial structure in the atmospheric 
signal. According to the Kolmogorov model of turbulence 
( Ihtarakii 1961), the power spectrum of the fluctuations in 2 
large 3 dimensional volume eh proportional to *?' 1 \ where 
q is: the wavenumber. The projected signal observed on the 
%ky cars have either a q 31 J dependence, if the wavelength 
is small enough that the turbulence can be treated an three 
dimensional, or a q * 1 dependence, if the wavelength is large 



DaL and Map Making 


5 


compared to the thickness of the atTnospheric layer supporting 
the [urbuicni motions {Church i 995; Lay & Halverson 2000D. 

Figun: 4 show* examples of the eunosphsric signature in 
"POD power spectra. Figure 4a displays the average TOP 
power 'tpccim from four groups of observations with the tele 
■scope not scanning ("stare” observations'). When scanring. 
the knee —where the power law metis the white noise level- 
inereasss in frequency by around l Hz and harmonics cf the 
scon frequency leak info the spectrum. Before averaging the 
power sjKCtra from different observations (data files), the av- 
erage power spectrum from the dark detectors was suterictEd 
to isolate the atmospheric signal from instrumental 1 // noise. 
TO Eh. fctre binned as a function of the FWV level measured 
by APZX during the period of data acquisition. It is ;1ear 
how (h= atmosphere signal grows with PWV, shifting ihe knee 
toward* higher frequencies. On the other hand, the power 
law ineex stay * rather constant near the 2-dime raw sal regime 
value The average pow er spectrum from the dark detectors is 
■■hewn for comparison. Ii is dominated by the thermal 9ue 
tuaiior . 3 of the cryostat. The fad that the “dark-dctccTor" 
power spectrum appears higher than the others in the fre- 
quency range between the knee and lOH/is ihe result of Ihe 
power subtraction mentioned above. This, indicates that the 
1 if pljjs rradoui noise dominates over the atmospheric plus 
detect it noise in that frequency range. The knee frequency 
ranges from I to 5 I Lz for stare observations depending oa the 
weather condi fion*. 

The departure from a pint power law shown in Fig. -la im- 
plies dial the power law index varies wiih frequency. By fit- 
ting pewer Laws in small frequency ranges, shown as dashed 
lines ir Fig. 4a, we were able to measure this dependency and 
group 'he observation!! as a function of their power spectrum 
slopes at frequencies near the knee. We found that, for similar 
PWV conditiortft, the slopes can vary significantly, with in- 
dices ranging between the 2-dimen&oria] and .1 -dimensional 
regimes, os shown in Figure 4b. Low frequencies which are 
related to large scales in the sky. are dominated by the 2- 
di meridional regime, with a power law index of appro xim sidy 
a = -2.4. Al frequencies around 1 Hz the power spectrum 
from some observation*- ibluc dots in Fig. 4b> follow a steeper 
power law. suggesting that, under curtain weather conditions, 
the 3 -dimensional re&iinc dominates at these smaller stales. 
Ihe slopes lend to increase with wind speed. Higher wind 
speeds are expected to shift larger feature* of the turbulent 
partem towards higher TOD frequencies. As large scales are 
domin/.ted by the 2-dLmcnKional regime, this would cause the 
opposite effect of reducing die slopes at higher frequencies, in 
contradiction with what \H observed. This result suggests that 
big he: winds might be associated with intrinsic properties of 
the turbulent layer, such as its width or height, producing a 
steeper power law (Church 1995; Lay ft Hal verson 2000). 

Atmospheric structure* larger (him 24' — fhe held of view 
Of a detector array- appear as a common mode among de- 
tectors, while smaller features can produce sub -army “corre- 
lated modes," as described in Appendix A. Given that the dis- 
tance to the turbulent layer is commonly less than a kilometer 
(Robson ct at. 2002), it becomes out of focus, smearing out 
to scales of nearly I0\ in agreement with our optical si mu 
(□Lions. This correspemds to roughly a third of the array size. 
Then, given the scan speed, such signals appear at frequencies, 
between 2 or 31 lz in the TOD. In general, the common mode 
of the .ielccuirs is a jumd estimator of the atmospheric sig- 
nal, but the estimate enn in principle be improved by dividing 
the army to account for sub-array ainm sphe rc s [juncture. Our 




t-'io 4, Armo-spheric ^ignaturr in I'OD power spectrum. (s) Av 
erage power spectra for four groups; uf non-&tanning TODs .>tdiacd 
by PWV level. The evenige power spectrum from iht dark deteefeffl 
was, subtracted from tsaeh gimip pnwer spectrum and Is shown for 
comparison. The dashed linesi shew a power Saw fit to each sjwh:- 
Lrum. The legend indicates the mean PWV and the power-law index 
from the fit for each group- 1 .ogaTitbmic tinning was u&ed to reduce 
i he noise in the plot. (b> Power law index as & function, of frequency 
Fir ihe average power spectra frysm three groups (blue (•), green (?) 
and red (A» Ls-f observations sclccied by dies: power law iodex at 
frequencies dune m ihe knee, all belonging to die diird PWV se- 
lected group in. (a). Each value was obtained by fitting a power Law 
to Uic corresponding power spectrum in a small frequency range (in 
logarithmic space); she. l isted frequencies are the mean valutsior ibe 
associated frequency range. The error bars show the dispersion of 
dtc indexes from die members rif each group. The same is shown in 
brown {•) for the- dark detector average power spectrum for compar- 
llKWI. 
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TABLED TABLE 4 


MeUlrJJ VALUE OP TIIS DTPlCAL atHSII IVI'rt I'Ll It OKlILCTOfc Ah It NOISE COUTH IUUTION SUMMARY. 

Til G TOTAL SF BS Itl V ITV HJH Ti l UAH HA V AT TWO 19 HOUKNCY 


RAMCliff AFTFR HT WOVISn fill' 3S MlHlK.S W LtJURI LATIN 

NDISli. IMF ! TN 1 L S : i U v,l I 1 l-MI’IUI ATI F. FOft CWtl 

N.ilhc Kourca 

*i\LV | mK* sE 


1 . 1 cruA'msN. 


tn-bim-0 (fe^cUnir noise 
Aliened noise 

030±(l. 10 
0 i7±o.ue 

Ha I 1 ?* 

NI Ti >t lf.K^i| 

NET^i liuKv'il 

Optical lollin': 

0,13 

s-atmit* 

100 - 120 R 7 

7hn I 2H 
1051 ±21 

317*3,7 
43-7 ±3.3 

Tbia) ndiw 

t>-«±0.04 


' Thc^ptl tmnd in < 30 K? no th^ qnlria hatcVllTMt* lh* Initntmeal 
scnssdvily 


otic nip's 10 detect coherent motions or atmospheric Icaiurcs 
across the array, as would appear in a moving frozen-sheet 
model of the turbulent layer were not successful. Instead, we 
suppress the atmospheric noise hy vilorii? for ilic strongest 
atmospheric modes in the time stream during the mapmuking 
process, as discussed in M L 

S. SYSTEM SENSITIVITY: UNCOR RELATED NOISE 

Above the almuspfiere knee frequency, the TOD j> dorm- 
Ttalcd by broadband random nos so this noise is generated at 
the detectors and readout circuitry and it is essential ly linear - 
related among tteteclorv After de- projecting the correlated 
mocks from the atrnusplH-rc and systematks (tee Section 6 
and Appendix A). this noise cun be measured by averaging 
the TOD power spectral density (PSD) over the desired nmgt 
of frequencies 

The total sensitivity of the array, expressed ait the noise 
equivalent temperature (NET), is given hy 


NETirt =• 



1 

NET? 

" J 1 


-i a 


( 4 ) 


where NET; is the NET of cadi working dclector. Thus, the 
typical sensitivity per detector can he dc lined as NRT,^, = 
NET W xA'it. where is the number of the "effective" de 
tetters, as defined in § 10, 

Table 3 lists the typical values of the total and average NET 
for a tnid- frequency range (5-3GHzj and a high-frequency 
range (10G-120HZ), in equivalent temperature units for CMB 
fluctuations, Unccnaintic* show the dispersion among 15 
minute TO IV The noise increase at higher frequencies ls 
driven by intrinsic properties of the bolometers, as described 
in Marriage (20061, Zhao cl ah (7008) and Nicinack ct ill. 
(ZOOS), However, our si £ 110 1 band is Siinilcd Ed below 30 Hi by 
the telescope scan speed and the beam size (see Figure 3), 
Thus, the mid- frequency range is the best e.stimnie of the in- 
strument sensitivity, 

As given in Table 4 t the typical noise variance per detector 
is 0.62 i 0.04 mKV. The total noise can primarily be sepa- 
rated into in-balld detector noise, aliused delect or and readout 
noise, and phuiun noise, all of which add in quadrnture. 

To reduce tile effects of noise aliasing, the deteeEorx are brst 
sampled at 15, 15 kHz, then a digital four- pole Butterwnrth 
low-pass fiber with a cutoff frequency of 122 Hz is applied, 
and finally the data are resampled ai 398.72 Hz. The detector 
noise bawl width is limited to 8 kHz by a 700 nH inductor in 
series with each TES. 

To determine the effects of optied loading and aliasing, wc 
performed a dark iesl in which wc opened the cryostat, put ;s 


4 K t reflective) cover over the detector 1 !, ^nd collected data at 
n variety uf sampling rales. The typical noise level in the dark 
was O.flfiiO.HTtiK 1 *. where the uncertainty shows the dis- 
persion among detectors. The unit?, are equivalent to the ones 
itl Table 3. After fitting die noise as a function of sampling fre- 
quency. die in band detector noise yielded 0 30 ±0 lOmK^s 
and the total aliased noi se contribution was 0, 1 7 ± 0.08 mK 1 s. 
The latter includes aliasing from both detector and readout 
noise. 

The readout noiic is dominated by SQUID noise and 
preamplifier nuisc. EhiUy sampled, the readout noise is esti- 
mated to be around 1 .2 x 1 0 1 mK 1 s ' based on 50 MHz ntca 
suremctiEs), which js expected to increase by roughly a factor 
400 when sampled ac 1 5 kHz, reaching Q.GSmK 1 ®. The 
SQUID noise aliasing was significantly reduced in season 
2010 by reducing the readout bandwidth. Taking all this 
into consideration , wc estimate that slightly more than half of 
rise total aliased noise contribution is aliased detector noise, 
which is comistcnl with the detector noise aliasing analysis 
presented in Niemack (2068), The other half is SQUID noise. 

The dark tests also revealed that the optical loading con 
tribute* 0.I5mK'\ of photon noise, By measuring the satn 
ration power of the detectors at different atmospheric condi- 
tions ihrougbo-ui the se^orq and comparing I hem to the satu- 
ration power in the dark, we found chat the optical loading is 
2-24 pW when the PWV is 0 mm. This loading is dominated 
by spillover contributions, emission from optics and dry atmo- 
sphere emission. The spillover contribution was reduced by 
0.36 pW in season 201(1 by adding a ba filing structure around 
I he sevondiiiy minor. In the same way, we determined that 
water in the atmosphere contributes another 0,7 pW/mm of 
loading, so in nominal conditions (FWV = 0,49 mm) the total 
optical loading i.*i 2.59 pW. This contribute?! to the noise as 
photon noise: fora fully incoherent detector coupling to both 
puLarizaihnis, the noise in units of power squared per ur.it fre- 
quency in given by 

NEF 1 «s 2kvP+-?z—, (5) 

where h is Planck’s constant, u = 149,2 ±3.5 GHz is the 
oahlral frequency 11 , Av - I8.4GH7 is She band w tilth, and 
P = 2.59 pW is the power absorbed hy the detector (Zmuidzi 
mis 2063), The first itrni in the equation corresponds to pho- 
ton shot noise while the second term corresponds tu noise 
from photons arriving in bunches, as is expected for thermal 
radiation when ihe occupation number is high. This leads to 
a photon NEF of 3,0 x IO~ l7 W/yi1z, or roughly O.J7mK 3 s 
in CMB tcmperulure units, in agreement with our measure 
metUft, 

11 The eencfttl frequency fur * kiiylrighdMni m-jnee. 
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Fin. 5.- Total array NET curing she 2tKJK wesson in the mid- 
frequency range (3-2flHz) for 148 GY\z, The values are grouped 
Li cays anc Uic error bars are the standard deviations within cadi 
day. The f’WV is al&D plotted. for reference (dashed line). Values 
are calibrated to CMB equivalent stniEs. Before computing ihc sensi- 
tivity, eight mulU ixmimon modes, four row-oomlsned modes, four 
column correlated modes, and the residtKi twelve modes with high 
cst ^rctgulSF varies were removed (see Appendix A for details). The 
noise improvement after September 24 came from mitring off some 
oscillating dcsocEors, which were ocmlamiinaling neighbors. 


Table 4 shows a summary of the different uncorrelated 
noise contributions Dial determine the system sensitivity. 

Figure 5 shows the mean noise between 5 and 20 Hz for all 
data files So the 2008 season. This is the tola) noise which 
goes into the maps and can be reduced only by increasing the 
observation time. Noise estimates from the maps can be found 
in Marriage et a!, (20 E ] a) and Das et al. (20 E 3 b). 

6. SYSTEMATIC NOISE 

In addition to the atmospheric and random noises, the ce- 
lestial signal is contaminated by several systematic noise 
sources. 'Ole important ones are thcnnal drift of the cryo- 
stat, mechanical accelerations which couple both optically 
and thermally (by causing detector temperature oscillations), 
electromagnetic pickup and magnetic pickup. All of these 
effects cause, zero lag correlations among different detector 
TODs, or can be directly measured using the dark detec 
tors. Detectors are distributed in 32 “columns" and 33 “rows," 
where the last row contains the dark detectors. Each column 
shares die same timcmultipcK-ed SQUID readout circuit (de 
Korte el ai. 2003) within which each row is read simultane- 
ously, Thus, lltere is one dark detector per readout circuit, 
serving to assess systematic noise from it. Moreover, the high 
redundancy of live detectors across the array can also be used 
to assess systematic effects (see Appendix A), 

6.1. Thermal Drip 

'Hie current signals from the boloEncters arc amplified in 
100-SQUiD series arrays (SA) operated at 3 K (Swetz et a I, 
2011), Slow temperature changes in the SAs produce stow 
drifts in the TODs. During the lirsi 10 hours of the eight, 
the SA temperature drifts down by 250 mK, and rises up by 
150mK in the East 2 hours. In terms of equivalent sky temper- 


ature at 148 GHz, this corresponds to a drift of nearly 4 Kin 
signal. Eh frequency space, the drift imprints a I// signature 
on the data, which meets the detector noise level at a knee 
of nearly I Hz. Despite small differences in the responses of 
different SA amplifiers, most of this signal appears as a com- 
mon mode to all the detectors. As die coupling occurs at the 
readout circuit, this signal is also present in the dak detectors, 
simplifying its, Identification and eventual removal. 

In most cases this signal is subdominam to the atmosphere 
signal, but when the PW V is low enough the thermal drift 
becomes comparably significant. This is not evident in Figure 
4a since the examples of low-atittospber ic power gel averaged 
down when grouping many TODs, 

In contrast, drifts in detector temperature cause no measur- 
able effect because it is servo-controlled to within less than 
1 mK error. 

6-2, Mechanical Accelerations 

Mechanical accelerations, occurring mainly at die scan 
turnaround, can couple optical Ey and thermally to die detector 
and result in an undefined instrumental response. 77ie opti- 
cal coupling can be produced by a mechanical motion of the 
detector coupling layer, which is a 40/im thick silicon layer 
placed lOO^ffi away from the detectors for optica] impedance 
mulching. The coupling efficiency is sensitive i_o the distance 
between the coupling layer and the detectors, so small vihra 
rio.ns can cause significant effects, especially at higher fre- 
quency bands where the coupling layer i a thinner and the cou- 
pling efficiency is more sensitive to changes in the distance. 
The effect is also larger near the center of the detector array, 
presumably because the coupling layer vibrates in its funda- 
mental mode. En the TOD. vibrations manifest themselves 
as a series of spikes visible in the common mode at the scan 
turnarounds, with opposite signs for opposite directions, lead- 
ing to tinea in the power spectra at several harmonics of The 
scan frequency. In the 14$ GHz “waterfall plot" of Figure 
6, the vibrations contribute to the scan harmonics and .some 
resonant lines at higher frequencies. The latter are also seen 
for stare observations, and iheir central frequencies are shared 
among all three arrays, suggesting that they correspond to nat- 
ural frequencies of the whole system. To mitigate the mechan- 
ical effect, the turnaround acceleration was reduced from its 
value in the 2007 season to the value shown in Table I af- 
ter October 8, 2008. The coupling layers for 2ESGHz and 
277 GHz were removed after the 2008 suason. 

Thermal perturbations of the cryostat. revealed In spcctraS 
analysis of the bath temperature, also add to the low frequency 
scan harmonies seen in Figure 6. This effect is column- 
dependent and differs from the coupling layer effect in that 
the signal is not stronger at the center of the array. 

In general, high frequency spectral lines in the TOD from 
different detectors destructively interfere when projected into 
map space, mostly canceling out when many observations are 
combined. On the other hand, low frequency harmonics of the 
scan may produce non- negligible bar- like features in fhe maps 
perpendicular to the scan direction, if not treated properly, 

6,3. Electro m a g netic Pickup 

Electromagnetic pickup couples to the readout circuit, pro 
ducing various signatures in the data. These signals ap 
pear correlated among subgroups of detector TODs, par- 
ticularly among detectors from the same column or row. 
We call these "column" or ‘'row" correlations respectively. 
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U'i BeEcur icmD'vuig cotreliied model t^f A fit/ ntmoviflg wjitc! iled rnKitv. 

Ftc 6 Frequency sp^-e wait: fait pkH cmiLumng the power ipttlruni frum all live detectors ffmm the HtiC-Hz band, before and after 
remnvrr.E 20 cemeigicd modes. This TOD was obcairaai under aomtliiJ '.cunning eomJiikitii (/**„ ■■ 0.1 Hr). on November 20. 2008. The date 
.are calibrated Ul CMB temperature equivalent ufutv bach huncMHii! line in line 213 plot represent!; the power fpoLinim uf a single detector. The 
black horizontal Lines separate detectors from different columns Mntl i>f the low !rcquaiey lines ure ir-un harmonies. They arc explained by a 
combinsiinn of effect* including spatial variations in the atmosphere brightness through the neurly triangular scan pattern, thermal n&cilktions 
of the cryostat amf coupling layer n.wHIIaiinsiK. Die thick lines around 5 Hr (soc columns 2- 5) are also present in start obsietvwisjnj, They are 
thought to be resonant frequencteii uf the system. Note (hat the science band i« between about I a nd 20 Hz. After removing 20 correlated modes 
from every detector TOD* harmonic tenures arc significantly reduced, ns i* the low frequency power fmm (he atmosphere arnl from thermal 
drifts. When making maps, the correlated modes are identified und properly tie ‘Weighted (o produce an unbiused map solution,. 


each corrc-sponding to a different source of electromagnetic 
pickup. Narrow-band signals appear correlated onion i/. Ji in: 
tors within the .'faint column. They couple in somewhere af 
re r il>e first SQUID stage of llic readout circuit, its the nub 
sequent circuitry is shared by the detectors in a column. Du 
the other hand, broad-band signals appear correlated umong 
detectors within the same row, or even a few rows apart. 
We believe these are caused by rapidly varying signals, Auch 
els spikes from strong current transients: given our time- 
multiplexed readout scheme, in which the same row is rend 
from nil cahjums simultaneously, these signal* become cone 
lated among detectors in the same row. 

figure 7 shows the detector correlation matrix for u given 
15-minute TOD fik, where every dement in the matrix corre- 
sponds to the correlation between the TODs of detec tors i find 
j, which index the dementi of each axis. For instance, the 
diagonal elements correspond to i — j, so they are all equal 
to one. In Figure 7a, adjacent detectors on each axis belong 
to the same column m the array, with a black line separat- 
ing detectors from different columns, while Figure 7b k or- 
ganized by row's. We can clearly see the kinds of correlauons 
described above affecting either columns or row* of detectors 
As would be expected from a stochastic distribution of tran- 
sients, the patterns seen in the row-correliited matrix change 


in H me. To quantify the correlations, a quality factor is dc 
fined as ibc rotten of I he squEircd off diagonal elements of the 
CQirclution atmtrix: 


fi- 


2 

A(iV |) 



( 6 ) 


where c, t is the correlation between the TODs of detectors i 
and j t and N is the number of detector TODs. Ijawcr Q ["acton, 
indicate less correlaicd noise. 

Ik cnlinnii rorrdnied signals were significantly reduced 
by cure fully shielding the system from electromagnetic noise 
in the environment. The row -cor re Lated signals* instead, were 
found to be reJmcd to the switching power supplies feeding 
the readout electronics. They were replaced by linear power 
supplies at the end of the 200$ season. 

If net treated properly, broadband row -correlated signals 
can produce features in the maps, mostly seen as features per- 
pendicular to the scan direction for that particular observation. 


6,4. Magturic Ptdtup 

The readout circuit of the detectors uses SQL 1 IPs, which 
arc very senxirive to magnetic fields, Despite significant mag- 
netic shielding, some residual pickup remains as ihe telescope 
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1c) Ctrijinn (i-rnmin:. Mifitotk pit l. jp, itmutplicft and dark muila nanuwd. Irt) How dorainant Magnetic pickup. atfTWtpbere and duk mods removed. 

14c. 7. — CutTeluridn TTiEtnots for a 1 5 minute dam stream ttaiirvcd Ifceemhtr ), ?0(]7 The 55 1 active detectors arc each nirmherea, starting 
from the upper left of tbc crime In: inn main*., m order of columns (left pitch) or mwi (right p*wl^. The remaining 47? Jriecrswri wot cat 
Iran this particular data stream duo to Inadequate performance (see See. 10.3). The solid linen divide the 32 individLial columns (left panels) or 
rows (right panels]; (he diffcrc nt widths between the wild line.* a- fleet the varying number of active detetinn in each column or now fcf. Fig. 8). 
The magnetic pickup and the atmospheric signal hove been removed from the time stream try subtracting a betf-fil sinusoid for the former and 
by removing eight mutt: -com men model for the Lsuer (see Appendix A). The top panels, show these (itters only, w uh a quality factor Q = 0,Q& I 
(Tq. 6> The bottom panels addi tUmidly have dark correlated imulr* removed, with an ohvknm Imge suppression of detector carrelaiian and a 
quality factor of Q = 0,015. 


sweeps through the Barth's magnetic field. Although the scar 
is almost triangular, live observed signal is chiefly sinusoidal, 
which can be explained by eddy current teases and hystero 
sis in the magnetic shielding. The magnetic signal amplitude 
and phase arc fairly stable across defectors in the same col- 
umn, but differ significantly between columns We believe 
these phase shifts are related to the complexity of SQUID* 
and electrical loops in the system. A typical amplitude of this 
signal is equivalent to nearly T —7 mK in CMB units. 

Ah this is a readout' -related signal, it can he measured using 


the dark detectors and removed by lilting and subtracting a 
sinusoid, Moreover, die scan frequency is significantly below 
the lienee band and the magnetic pickup is heavily down- 
weighted in the mapping. 

7. CALI EHRATIGN 

Ottr calibration takes inio account the detedur properties 
and electronics; camera and telescope optics; and atmospheric 
conditions. It consists of three steps. The lirel adjusts fur vmi- 
alkmji m (Selector rcfiponsivity due to changing atmospheric 
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luaduip iinJ aimo.vphcnc opacity. Next, a dctecror-ckitcciDr 
j^ain ie*tn is L ip plied la account for relative variations in opti 
cal erupting *ia:us> the camera. Finally, j single iKHTn.-I.iza- 
tkm is applied lo all the data to account for the overall sir item 
efficiency for celestial source* This procedure is si.milir to 
the ant presented in Swu/cr (2008 J. 

7, 1 . Calibration Vanattom vrith Time 

The nightly variations in aimosphcnc loading and daily 
cryogenic cycling necessitate a ic biasing of the deteetcis at 
[he beginning of each night's observations, which affects the 
detectrr response. Additionally, the opacity of the atmosphere 
varies from nighi to night. Changes m loading and opacity 
throughout til* night produce smaller variations in the system 
response, 

In order to account for the variations in system response 
between night*;, a rraponsivfcy calibration is derived from the 
analysis of 1-V curves taken during the nightly rebi using. An 
/ V curve is the relation between I he output current response 
(/) of Lhc detector's .SQUID readout circuitry to a slowly 
tamping detector bias voltage (V). Por the 148 GHz band, the 
nlrdiiLr'i deviation in rcsporisivity between nights as measured 
by the 1 - V curves is 3.0%, 

A second test, known ns a bias step {Nit-mack 2008), is per- 
formed three times per nig, lit to detect changes in detector re- 
sporisivity through (tic night. This is achieved hy recording 
live detector response to a series of small, square-wave pulses 
applied so the detector bias voltage. The responHivities can 
be obtained by analyzing Ihe detector responses to ttiis signal, 
and nre found to be in agreement with respomivities obtained 
from the i - V curve analysis. F'orlhc 1 48 GHz band, she me- 
dian deviation in icsponsivity over u night, as measured by the 
bins step, is 1 .0%. For more information about the ACT detec 
tors* the biasing routine and responsivtiy, sec Fisher {2009),, 
Huttisidli et id, {2008), Swell el al. {201 1) and Zhao cl al. 
(2008) The time dependence on the atmospheric conditions 
will be discussed below mgelher with the overall system cali- 
bration 

7.2. Time Independent Relative Detector Calibration 

To determine the detectors’ relative gain cofFliciunts, wc 
use the large common-mode signal provided by Hie variations 
in aimcsplieric brightness This is analogous in a flat-field cal- 
ibration done with an optical CCD. For each 1 5- minute datu 
tile, wc compute the gain factor that best fits the detector drift 
to Ihe common mode drift The dispersion of this facto: per 
defe c irr over the season, averaged over all detectors, is better 
than 2% rms for the 148 OH/ array, This includes the variabil- 
ity of the time-dependent calibration step. Wc correct fur this 
by multiplying each detector time stream by its corresponding 
gain fa dor. 

7.3. Oventit System Calibration 

An rvernll system calibration is established by rmirf 
merits of Umnuv Hy comparing the peak response of the 
instrument to the known brightness temperature of Uranus 
(taking into account Llw ratio of Uranus 1 solid angle, fi L = 
0.235 ~ 0 OlOnsr. to the telescope's solid angle at 148 GHz. 
ftjHCT - 215.8 ± l.5mr). an overall calibration of the system 
response is obtained. The statistics] error derived from 30 
measurements in the 148 GHz hand is 1%, while solid angle 
uncertainties are below 1%. The equivalent Rayleigh- Jeans 
temper it Li re of Uranus is obtained by reprocessing the data 


from Griffin & ckton (I993) t in combination with WMAF-7 
measurements, of the bnghtnes* of Mars end U ranus f Weiland 
etaL 201 1). Tt yields Kft,5±6K al 1 4£ GHz, producing a 6% 
uncertainly in the overall calibration. 

The conversion to sky temperature al&o depends on the at- 
mospheric transmission. ITiis is estimated with the ATM 
model iPaidc ct al. ^10!), which uses PWV measurements 
from APEX corrected for trie ACT sire and os her Atsca ma- 
specifk parameters. When PWV measurements are not avail- 
able, the season -average transmission (T = 0.976) is used. 
For ihe 148 (9fa band, the rms of the transmission during the 
2008 .season was 2.3%. Given that large atmosphere loading 
can excite non-ljncarilie£ in the detector rrspons: vines, an cx 
tra degree of freedom is given to the planet ill as a function of 
PWV. producing a final correction to the calibration factor. 

The overall map calibration is compared to the WMAP-7 
year map, by correlating mutiipofes in the range 400 < £ < 
1000, providing an uncertainty of 2% in temperature (Hajian 
el al. 201 1 j. This analysis shows that the Uranus calibration 
is a factor of 6% lower than the WMAP calibration, which is 
within ihr expected uncertainty. 

Purring 'ill pieces together, the overall system calibration 
for 148 till/ in the 2008 season is 

F{w ) = C ex p [fo + T* hm- t x (w - w)) / sin „ (7) 

where C - 19-41 K/pW is the overall calibration factor, fiU 
is the observation altitude, r rf ^ 0,0093 is the "dry opacity," 
t w = 0.0190mm 1 is the "wet opacity," = 0,0138 mm -1 is 
the the extra degree of freedom, w Is the PWV measurement 
from APEX, corrected for the ACT site in millimeters and 
w - 0,44mm is a pivot PWV used in the fit, Here Tj and r w 
are fixed parameters obtained from the ATM model, w r hilc C 
and Tj ore free parameters. If r„ is allowed to vary, we recover 
the same value That is predicted by the ATM model (Lhncks 
2009). 

8 POINTING 

The pointing solution is decomposed into relative pointing 
between detectors and absolute boresight pointing, both of 
which we discuss here. 

8. 1 Relative Detector Pointing 

The relative pointing between detectors was determined by 
modeling the beam in the time streams of planet observations. 
This analysis made use of approximately 30 observations of 
Saturn at the normal CMfl observing altitude of 50? 5, in Iwu 
a/imuib ranges corresponding to the rising and setting of the 
planet. 

The telescope scons are rapid enough that each detector 
samples ihe vicinity of the peak response io a planet several 
times in a single observation. Ii is thus possible to model the 
beam position and shape in two spatial dimensions, project 
thus to n detector signal using the telescope encoder informa- 
tion, and fit the data in Ihe time domain. For simplicity. w* 
use a 2 -dimensional Gaussian as the beam model. Tbc fit pro- 
duces azimuth and altitude oil sets fur each detector relative 
to the telescope pointing encoders, as well as measures of the 
peak response,, beam FWHM, and optical time constants. 

When combining different planet observations, wc fiat 
aligned [hem by applying an offset correction to each one. 
The position of each detector is taken as the mean of the posi- 
tions from all observations, after rejecting outliers, This pro- 
duct* relative detector positions for each array that are in good 
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agreement with design expectations (Fowler er al, 2007), A 
comparison between the offsets for rising and getting observa- 
tions {'vbich differ in azimuth by approximately 95 D ) shows 
no significant rotation nr shear of the array relative to the 3o' 
cal altitude and azimuth axes. This constrains the till in the 
telescope azimuth axis to be smaller than V and indicates that 
rotation of the army with respect to azimuth is negligible. 

The uncertainty in the relative pointing is no greater than 
1 "5 for all three arrays. Given ihe large number of detectors 
in each array, this small uncertainty in the relative deiecEor 
positions amounts to a negligible contribution to the pointing 
error ir. the final maps. 

S.2. Absolute Boresight Pointing 

The azimuth and altitude of the telescope encoder orbing* 
must to corrected to account for their offset from the true 
boresight for each frequency band. The correction is differ 
cm fur each of the four central CMB observation azimuths, 
namely the equatorial and southern stripes ai rising and set- 
ting orientations, 

The correction is obtained by projecting the data from t par- 
ticular band and orientation Into a map (after having applied 
the relative pointing solution described above), and compar- 
ing the positions of the bright point sources to catalogue posi- 
tions, The resulting offsets in equatorial coordinates are con- 
verted ic offsets in the boresight position, and this procedure 
is iterated to ensure convergence. The boresight offset varies 
by about I 'over the season, primarily in elevation. 

The dominant .source of error in the offset correction conies 
from estimating the centroid positions nf the point sources m 
die maps, before matching them to the catalog positions. This 
error scales inversely with the square root of the number of 
point sourees available, which was 20 for the 1 48 GHz band, 
and with their .signal Eg noise ratio. The uncertainty was found 
to be 2. M 6 for the southern stripe and 5 "3 for the equatorial 
stripe at 34S GHz, These values were obtained by adding, in 
quadrature t he error in the fit from the rising and setting maps, 

The telescope pointing is expected to vary slightly due to 
thermal deformation of the mirror structure. This variation 
is estimated from observations of Saturn taken at nearly die 
same azimuth and altitude on different nights. The rms poinl- 
ing variation of such observations is 4"3, Since any pixel in 
the final season maps contains contributions from many dif- 
ferent nights, this random variation docs- not contribute signif- 
icantly to the pointing uncertainty. 

Pointing deviations are significantly higher at dawn when 
the telescope temperature begins to change more rapidly, 
showing a trend that repeats every morning. Hie drift begins 
nearly } hour after sunrise, reaching nearly 50" an hour later. 
Data teken more than I hour after sunrise are not included m 
the maps. 

On top of this, the attitude was observed to drift by about 
20" over the course of the 1 20 day season, producing pointing 
errors that were correlated with right ascension, To remove 
this trend, linear corrections of 0"2/day and 0," ] 5/day were 
applied to the rising and sorting fields respectively. These cor 
rcctions have hide effect cm the maps other than to remove the 
small, RA-dependent pointing offset. 

The net pointing uncertainly is thus dominated by the sys 
tematic uncertainty in the alignment of the rising and setting 
maps, ^long with some residual variation due to temperature- 
dependent mirror deformation. Comparing the positions of 
bright point sources in the maps to their catalog positions, 
wc estimate the pointing error in the final maps to be IT 8, 


with no preferred orientation (Marriage ct al, 201 Ja). Mote 
that this uncertainty' is much smaller than the beam size, and 
would cause errors of less than 1 % in the measured Dux for 
point sources. 

9, DETECTOR 71 ME CONSTANT 'S 

A detector's time response is limited by its electro thermal 
properties. We model theft optical step response as an expo- 
nential decay with a time constant t , A finite response time 
results in a small shift in the spatial position of a point source. 
The shift depends on the scan direction. In the fits for the rel- 
ative pointing solution described above, time constants were 
i fit Ended in the beam model a.% a single -pole low-pass filter 
in the time domain. The time constants arc measured with 
a precision ol'^T < 0.5 ms, which is an upper bound to their 
dispersion from the analysis of 30 Saturn observations. Tint 
modi io time constant of the MB GHz array is 1.9 ±0,2 ms 
(f&g = 83.8 ± 8.8 1 l?X with only a handful of detectors show- 
ing responses slower then r = 10 ms {fat y ss 15 Hz). For com- 
parison, Equation 3 implies that £ ^ JO 1 corresponds to 27 Hz. 

K>. DATA SELECTION 

Data selection can be divided into two types: date file selec- 
tion (all detectors) and single-detector TOD cuts. The former 
determines the number of observing hours, while the latter dc 
[ermines the number of “effective detectors'" within each data 
file, defined as the sum of the fraction of die time that each 
detector passed the data selection: 

i?K4 

% = ( 8 ) 

7 ^ Foul 

where J] us*;* is the available time for detector i after data se- 
lection, and is the total time before data selection. The 
number of array-wide effective detectors can then be defined 
as the average of W eff uver all the available data files in the 
season. 

The remainrlcr of this section presents the data selection 
methods and Statistics for these two types of data cuts, 

10. 1 , Detector Classification 

The detectors are classified in three groups: 

» Live Detector Candidates: Those that have the potential 
to be used for map making;. 

» Dark Detector Candidates: Those that do not couple to 
the sky (for example a broken pixel) but can be used to 
diagnose systematic^. This includes the subset of the 
32 dark detectors which work properly and the subset 
of defective detectors with working readout circuit; 

* Broken Detectors: Those defective detectors; that can- 
not be used to probe systematic effects. They include 
those former live and dark detectors with, defective 
readout cireuits* and slow defectors (r < If) ms}. Many 
of them must be turned off while observing to prevent 
them from interfering with other detectors. 

Several metfteds were used to classify detectors into these 
categories. These include assessing correlation with the aimo' 
spheric signal {by far the hugest signal), searching for consis- 
tently oscillating detectors, and finding biasing problems. 

- 4 Wc have implertienred juUomatic methods lo idcEiSity n-sciSSating dc:ci> 
5SMS on sEtc Sty asd dlscomect them. 
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Thblo 5 give* 4 gummuty of the number ol' detector* in each 
of the three grot) pH for each array. 

Once tile live und dork detector candidate s have been iden- 
tified, t number of possible pathological behaviors may mill 
juiuify removing pan of the data from the redaction pipeline* 

10,1 Data Ftie Selection 

Out of the total number erf data file* acquired, we rejected 
file* lor the following reasons and in the following order 

* I lln IM cucropond if] planet obscrvaiitjnii Jind txher 
calibre tmn or engineering tests: 

* Data taken more than one hour after sunrise. to avoid 
pointing and beam errors caused by telescope riefornwj 
iii ms as a was thermally idling; 

* Inlet with fewer than 400 effective detectors, ns they 
were considered likely to be pathological; 

* Bad weather: PWV greater than TO mm finuwninsiwi 
belnw 911%); 

« Poor cryogenic performance: Defector base tempera- 
pin our limn ?mK above the nominal temperature, or 
when it changed mom than I mK within the 15-minute 
fie; 

* Foot calibration: If die relative gain dnperuHi of the 
detect nr* was more than 10%; 

* Dai* for which the analysis software failed. 

The unj] amount of dua available for analysis is wmmu 
rircd In Table 6 
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10J* OeteaorCutM 

Detector* am affected by sporadic pathologic*. Depend] rig 
on the kind of pathology, it may be necessary to reject a sec- 
tion or the full length of a detector TOD from a given rfoin 
hie. 

The mum causes for these pathologies arc quantum jumps 
of i br magnetic dux of a SQUID in the readout circuit f X-& 
jumps); excesove detector noise; conducted timse from oscil- 
lating detector*; excessive eketro magnetic pickup: and me- 
chanical contain i nation, which can he optically or thermally 
coupled into the signal. 

Hue following tests were performed over 1 5 -minute data 
tiles to detect these pathologies: 

* Drift teat: Ibis probes low -frequency devotions of a 
deieciorTODfrofn the atmospheric signal, The data arc 
tirsi low pax* filtered above 50ml 1/ and calibrated into 
unit* of power a* described in #7. Then the thermal 
drift in removed by dc projecting both the dark detec- 
tor common mode and the housekeeping temperature 
of the detectors. Finally, the atmosphere signal is re- 
moved by dc-pmjectiiig eight multi common mode*, as 
explained in Appendix A The dnfl error is defined as 
Live standard deviation of the residual data. Detector* in 
1 'IS GHz an’ eui as outlier* ii (heir drill error is higher 
than 0.35 fW. 

* Correlation test: The Drift test is complemented by 
finding the cmrclunnn between the detector drift and 
umty common mode drift for TOD frequencies below 
5UmHz. Detectors that correlate less than 98% are ex- 
cluded. 

■ Gam rest: The Drift test assesses the shapes of the 
TOEh. bul not their amplitude* The relative detector 
gain can be quantified by the factor that bent fits the at- 
mospheric drift to every single TOD at frequencies be- 
in'* 50 mHz. For thi > the atmospheric drift i* estimated 
as die difay common mode alter removing the thermal 
coniamitudoa. Detector* Me eat whene v er (heir gains 
differ from the mean by more than 15%. 

■ Mid frequency noise test: Some pathologies, mainly 
associated with mechanical cofUannruLnon. become 
prominent, si frequencies between 03 and 1.0 Hz. be- 
low the “science hand To tsoLrc pathological oetec 
tots, we band pass filter the data below and above those 
frequencies, de- project the array vuinmcn mode (also 
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fltcrcd p. and obtain the rtjmitird deviation of (he tend 
mb. which we call mid- frequency enur, Deletion an 
cut whenever thetr mid- frequency error U inure than 2-5 
time* the mcdicn value for the dam file 

* ihyh frrtjiiem v noise rem Ai rretpjcncics above 
l Hi the duui start being dominated by detector wise, 
which lit chiefly Guimian. Non-Gmutkianitiea above 
this frequency are mostly caused by electro magnetic 
pickup, conducted clcctriail noise nr opto- mechanical 
purturbations To isolate them, we high pass flit nr the 
dm n below 5 H/ and then tesl for ncm-Gauaiii*niiiea by 
computing die standard deviation,, skcwncaa and kudo* 
tii. Thi* ii done within secikrna of the TOD of unc 
scan period, with a characteristic length of 10 2 a. Us- 
ing die inns ftrnnat tun proposed by D'Agoutnn 4 He- 
UntcrCI990j, the hut two itehstic* yield a norma; div 
n but ion in lire cuve of Gauwian noise, which is veri- 
fied for mewl of our data Outlier* are itknli h cti by how 
much they deviate from the mean Only section* of ihe 
TOO with noise rms lower than 0-9fW are kept 

* icon ffji Motjnti defects and encoder errors are ajso 
detected and cut. affecting sub- section* of the TOP 
This is done by CAsrmnmg the azimuth nmc strewn and 
"Core tun# for interrupt ion* m tire scan pattern 

* Glitch test: The data ore also affected by spike -Itke 
glitches, fur example from cosmic rays Spike* larger 
l'kim in nmes the muse mss Of* cut. Including a it 5l 
hdTcr on either side. Also, whenever two such cuts ore 
separated by lew thin 5 s (hey ore stitched together miu 
£ single cut If more that 5 gtiiches are found in a single 
Crtrctut I 0D, then the whole detector ‘TOD i» cut 

* CalibtmiOn Bolometer test : During 2008 obKCWHtiocu 

* calibration boftotneter »u used to load the detecton 
with radiation of nrjgMy 400 mK every 24 minutes, 
each event lasting for 13 Bar tbcac, i window 

cf needy 3 s is excised around the event ' 7 Wc included 
this within the 'Glitch’’ cuts in Tibi* 5- 

As a general rale, if more than 2Q5fc of the defector TOO 
would be cut. then the fun detector TOO is cut instead* 

Dart detectors were selected in an analogous but simplified, 
way. In this case, only full detector TOD cute war performed. 
The cu: criteria were the drift error of lire dark defectors, the 
gain with respect to die dark common mode M and the noise 
rms. ui: given in raw date imiiis. 

If) 4. Data Sc tertian A’nwh* 

After applying nil (he selection enteria given mlinve. and 
considering only the date files available for analysis, the as 
erage number of efhetho detectors was 593 i 95 in the 
feawn at 1 4$ GHz The error here is the standard deviation 
over [jrffcreni date files 

Table 5 shows a summary with the number oi detectors in 
each or the three detector groups for the 14MGH* array, it 
well a* the cut contribution from die eight main selection cn 
tem: drift gam. mid-frequency noise, noise nrn, n vise skew 
ness. n>ise kuncwis, and partial cuts with more than 2fTfc of 

*’ the I'fJiSnjimn tatomelrr wn mu L-.-ii utiiir 2CXW 

** Sivicc Ihni Ihf ammtm mode or tbc dark dcLeGLitn la Juhti i-y me 
•kfnnil drift of the cryniW, which il the second lirjnl ii|iul in Uk live 
4c icc tor Jain. 


1 4# GH* uny. trtLi.l -.| »rerc| x>cim a jingle CetccSof. They 
■re oriented oi die may is prajeclnl on the sky TJotice dial some 
rows uni colmniaiffl: omtyseut. which IsoMUy dneto problem* 
in the tuning (rows) and readout circuus fcolnmnii. 


He* 0 ■ EBwtivc number of deiociori in Uk 148 GHz array in the 
2008 season. Circles denote duty averages and the error bars ere the 
standard! deviation wiihin iK*t day. The PWV level li ihcwr. with 
(ho dashed line. The increase tn the number of cfTerirvr detectors 

MtariLtig in September occurred tfter ■‘.■me ovdHailtig doecinfs, 

w hLch had betn eon tami n* ring other dcieciors. were tur ned olT. 


(he TOD cur Tlie bottom line is the average number nf ef- 
fective detectors m (he season l 3 igune 8 is a diigram of the 
14&GHz array showing life fruetton of Uk lime (hut each de^ 
lector u uncol Figure 9 ^bow s die daily tiumber of the effec- 
live defectors during the 2008 i»ca-son. along w ith ihe weather 
conditions indieated by the PWV. 

On top of these cuts. ■ imall fraction of the remaining data 
is cut during (he mupmiking process, as dcjicnbed m the fol- 
lowing section. 
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11 MAP MAKING 
II I Mapping Ettenlials 

In it nkc maps from ACT data, we solve for ihe best-fit sky 
given tic noi&c in the data. In particular, we find flic sky map 
that minimize* x 1 given a model l or ihe noise, and a node I 
for whul i Ik data should look like: 


dNMx-tn, (9) 

Here x in the model for which we wish to solve; M describes 
how th; data depend, on the model, and n is the particular re- 
alization of the noise in the ACT date, Trmiiiionally. i In a 
vector whose components urc the sky map pixels and M is the 
pointing matrix. In its simp Test conceivable form, each data 
point sre, 1 * a single pixr] in ihc map, so each row of M (eone- 
sponding to a single data point) has a single I in the column 
corresponding tolhc map pixel m which it was pointed, How- 
ever, out model for the data, x. can contain many contributions 
in addition 10 ihc vky inapt components include atmospheric 
noise, corral nied si gnu Is in the data, point source fluxes, and, 
missing (cut) data sample* as discussed below. The mapping 
formalism easily generalizes to cover multiple component* as 
[oug to Ihc data depend on them linearly: 


Mo[M,M a M, 


( 10 ) 

nn 


Wilh an estimate of the noise covariance M = (ua r ), which 

is a by matrix fo Ml mie - 10 s h we wish to find 

the model that maximizes ihe likelihood function 


£ = exp 


^d-M*) r *r ] (d-Mx) 


( 12 ) 


where d is a vector containing all the data sample*, The dan 
danl linear lea^i-squares solution is 


M r N“ l tV1x = M'N _i d. (13) 

The nufrix M r N _t M is loo big to be practicably inverted di 
recti y (we typically have 10 7 t nap pixels), so we instead itera- 
tively solve the least-squares equation for x u*mg a Ptecondi 
honed Conjugate Gradient (PCG) scheme (Wright ct «L ! W6; 
Hin&haw ct al. 2003; Press el al 2007), Preconditioning in- 
volves introducing P. an approximate inverse of DPft'M, in 
order to speed up convergence of c lassie Conjugate Gradienl. 
and solving the better-conditioned system. 

PM r N _l Mx^ PM r N f d. (Mi 


lb map ACT data, we use the mupmakmg code Ninkas t 
and run cm ihe Scinet Central Purpose Cluster (GPC) (Loken 
« at 2310). ITie mapmuking algonlhm for ihc maps in ibis 
release is improved from that used in previous papers (Fowler 
eta]. 2010; Marriage et al 201 la.b; Dunk ley cl al. 2011' Das 
ct al. 201 lb; Hajinn et al. 201 1 ; Sehgal e| al. 2011) in four 
ways; 1) rather than sol ving foi detector-correlated noise (in- 
dud! ng atmospheric noise) exp lie illy, we put the correlations 
in the rx>:se matrix; 2) wc explicitly solve for eul data (Lime- 
stream gap>); 3) wc 'iiibtract models for ihe point source 1 , m 
the ACT maps directly from detector rime streams; and 4) we 


re- estimate the nnist after creating an initial map to remove 
signal induced bias in the noise estimation. 

The mapping is done in a cylindrical equal-area projection 
with a standard latitude of f =-5 3 & -5 and pixels of 30" x 30", 
roughly one third of the beam FWHM 

1 1 . 2 , Pre-processing 

Trt addition to data selection, calibration and pointing, there 
are a few other pre-processing steps that are done to the date 
before solving for foe maps. 

We remove foe median of each detector time stream for 
each IT minute period, and subtract a single array-wide slope 
across die period ho foy* the ends of the time streams roughly 
Line up. Wc do this to reduce ringing iu Fourier transforms 
and to facilitate searching for correlations among the detec- 
tors. Wc h nearly interpolate across gaps in the data, such as 
those wiring from cosmic ray hits (again., to reduce Fourier 
artifacts: the data In ihe cuts are otherwise nut used). We de- 
convolve the time streams by the anti-alias filter (descrihcd in 
§5) and the detector time constants (discussed in §9). 

Next, the non -optical contatrii nation signals are reduced us 
mg foe dark detectors. This is done as follows: we take die 
time streams from each dark detector and subtract a mean and 
a slope. Wc then high pass filter the dark time streams. 1 ti- 
tering out any signal below 5 Hz, Most of the remaining sig- 
nal is in only a lew independent modes. Wc take those cor- 
responding to the seven largest eigenvalues of the resulting 
covariance matrix and do a linear least-squares fit to the live 
detector data, which had been similarly-processed (having re- 
moved the mean and slope, and high-pass filtered). Finally, 
we subtract the reconstructed fit from the data. We can do this 
because the modes subtracted are not correlated to the sky 
signal 

In Fowler et al. (2010), wc downsampled the timcsticamis 
from '100 Hz to 200 Hz, using a time-domain triangular ker- 
nel of the form [0.25 0.500-25]. The ACT si pal band (nearly 
30 Hz according to Sec 3) is well away from the downsam- 
pled frequency limit. Wc find that downsampling does af- 
fect the raw power spectrum at up to several percent for 
£ >~ 5000. I his should be mostly corrected for when us 
mg a beam measured using maps derived from downsampled 
data However, we also find that the SNR on point sources is 
2**- higher for non -downsampled maps. Consequently, wc do 
not downijnipk the timestreams tn the results presented here. 

Also in ccKitnist to Fowler et al, (2010), we carry oui one 
further step in Ihc time domain: With ihe 30" pixels somc- 
w hat undeJMmpling the ACT I .'37 beam, we find foat T to get 
percent level accuracy in source fluxes. wc must deal with 
the bulk of the source flux directly in ihc time streams. To 
do this, wc make first-p.is; maps in which saurce-s are found 
and their fluxes estimated, Simulations sfoow foal these fluxes 
are typically accurate to i few percent, with source flux sys- 
tematically underestimated by approximately 4%. Wt then 
take these source fluxes and model them directly in foe lime 
swans using the full ACT beam, and subtract the model front 
the lime streams. We do a simple source-only projection into 
a map. and mid this to the i mostly) source-free maximum like- 
lihood map Wc find that., in simulations, this recovers mean 
source luxes to I s accuracy, at which point it is suhdominaitf 
to the calibration uncertainty, We derno additional filtering in 
foe rime domain. 
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1 13 Naat Wodthnx 

TIk notM izrocture ts modeled in both Irequctxs and tirr< 
dunum* We include, as a term in the none, a lime d* tiwn 
windowing of the tirsi and Last 20 seconds of nxh Imie *3 ream 
of I he fuftti ( I - coin .i > / 2. again to rcdu ce 1 v kit irr migi n> ^ r 
Mien search for correlations acro&s the data by tin mi mug their 
enva rianee malnec*. split at 4 Hz U zr 1500). Wc find ill 
eipenvectors in the low-frequency covariance matrix with enr- 
responding eigenvalues greater than 3.5 J times the median 
(Iho eigenvalues of (he data covariance matrix A* A are the 
iqtunn of the corresponding data singular valuer)* to tlika CGT- 
respnud* to an amplitude of 3 5 times the median in the time 
smarm ) Vri* itteii project those eigen vectors out of the high 
frequency covanance matrix. and again find ati eigenvectors 
wiih eigenvalues. larger than 3.5 : rime* the median m tf e re 
irumng high ’frequency manii We find that the nMf« are 
mu pjdit ulirtv sensitive tc the exact AmfhOld v.tlucv ch t\c i 
Thti procedure typically firak 15 to 20 low frequency e.gen 
\Tcim and one or two high-frequency ones. The eifeavcc- 
ton are reaponc panerm of correlation* across the am v end 
ndy correspond id like a common mode, ptdkm 
etrn the array, and the row -correlated noiise. They provide 
ihe I incut ewnbiimtiOd of detector time stream* needed to pro- 
duce a correlated mode, represented by the vtciot i til I qua 
[ion A I fins projection corresponds to removing arour d 2t > 
mode* out of the Approximately 600 available. 

With ihe shapes of the army correlations In band, we can 
uac Micro In complete The description of the noise. Wc solve 
fur the correlated modes corresponding to the eigenvector* 
ei in I vl hi rni i them from the time streams. Wc then linrier 
mills In sin ihr- (dc-cunr luted) detector trim: -.i ic.un ni'J 1 1 w 
cone Luted modes, arid model them using Frequency biro. In 
each frequency bin* we hod the average power m each cciec 
tor lime xtream and in each correlated mode. If ihe detector 
nrim ore dewxesi by the diagonal main* N w , the sktcutuf 
covariance eigenvector by the ji a~w by jb^o. matrix V and 
their cnrre&pondijig tun-wise noise power* by Ny, then the 
bin wttc Hourrer- space noise is simply + VM V r 

Here N u in n«*vr by n^cu niZiti . x acting on ■ tragic 
freqi>r:us tun. Ibis matrix can be quickly roveried iring 
die Sherman Woodbury formula (Duncan (1 W). %« Hager 
11989) for a review). If we denote the time -domain edge ta- 
pering i ■ 1 ‘c-r ,Lt: u tis W ..r.J the Fcurrc: transform ■ ipre*!i n K 
F* then our entire noise inverse becomes 

NT 1 == W T F t Nf J FW, (1 51 

where the operator Nj- 1 is conformed by all the N^ 1 Much thill 
it arts i in every frequency bin in the Ixmritt irniUlfbrm In 
the pftviuiu equation, every operator can be represented by 
an by n^, ,,. i i.ieiix, in which detector operation* like 

\ f are exputidedxuch that all sample* within t single ds 
MOT are treated in the same way. To mtmiiii^ nm 
ttvity to iiny low frequency non Gausvan component cf the 
■tmtMpherlc none, we taper the Former weights at frequen- 
l'ics he'ow b ' Hh 1 , with the '*Cighi eiplszifly *-ci u> MB MP 
0.25 Hr We note that the noise in the mapciukin; L-tjuiii m 
(Bq, 14) can be interpreted as x set of weights The reap n up- 
UMtii i! tht weights are perfect, but thr map rcfiuin^ unf' jv ./ 
as long ax identical weights arc used for the left and ng hi ude* 
of the equation, and t he weights are anconelited with the iky 
signal Thi\ is the rvvrnliil pan of the meihod. Time ^Lirum 
fillers arr by design biased and must be aceouiilrtl tin with 
sim ul m ions. On the other hand, in our ireamnent, modes in the 


nup that might cause prublema are de- weighted as opposed in 
being hhofed, producing un unbiased wkiiinn through careful 
ronidcratkii of (he noise qroclure of the data. 

One tuTther difference remain* between these maps and 
those used m Fowler et li (2010) Mete, fnr every single time- 
xiieara sample cm, wc fit for it* i unknown) am pi nude ad pan 
of the map solution, «> suggested by s fur example, Pittanchon 
ct at. (200ft), In adduion lo cub from cosmic rays and tlie 
like, we explicitly col ihc fil’d and In , hi secund of each TOD 
to give the mapping algorithm the freedom to match the TOD 
beginning and end ns smoothly a* |«,>ssible ,Since thexc dam 
are already highly downweighted by the time-domain taper, 
the additional amount of dnLi lujts is ne eligible 

1 1 ,4. Map Solution 

With these pre processing uepv form of the noise, and 
form of the solution specified* »r muu ihc!i j^iually sohe 
for the map. We use the dank Frcconduinnrd Conjugate 
Gradient algorithm, using the hi I count map ■■ ■ Jacobi ! di- 
agonal) precoBiltioiKr fee the iky map pmt of the cokobaiL 
To recover the scales around t ~ 2)Kl and Ixlow, w t need a 
few hundred PCG imraitinm. The map mdueg; is quite CPT ■ 
i ntepuj vr i*ii±i each itentinn inking about a rmoute on 1600 
2.53GHz Nchflkm cores on ihc Sane) UPC, t>r a full run !□ 
ItKiO iterations teltinp atHHU n diiy nf wall dock lime, or 5 
year* of CPU time. 

Wc note that care must be taken when csiiiujimg the noise 
to make, sure it does not bias the mnp. Cun rider the follow- 
ing thought experiment: Two dnicdorN are sennnod across 
u source, and their weights arc civil mated from the internal 
Hcaitem uf their time streams In gnnentl, nolvc in (be de- 
tet'iivs will lead to one of them observing a lower II ux from 
the ‘iource. If the detector weight* art then measured am 1 
the twite ts not removed from the tunc %i rearm, the detec- 
tor which happened to the lower flux will uwafiy 

have a lower vsimcc in its data* and will therefore on aver* 
age receive more weight than the otter defector When the 
detector* are combined* the more heavily weighted Lew-dux 
rarovuretnent will lead to the voucce flux bang xyuemuicalty 
biased low. The magnitude of the bin n set by the SNR 
in the indmiM time vtream*. tM the final SNR The bias 
can he mitigated if an estimate of (he signal u removed from 
liic time streams before noise cftiiraitioii. Since the UMB i* 
highly subdominant to the imixe in ACT Umc riieam*. this 
efTcct ia amall. but simulation* rimw tiini i) is ru-r negligible. 
So, ftrsE we make an initial mup using thr full daiiivti tutd the 
noise measured directly from (he time tirumi*, iheii we mb- 
trucT that mnp hwn rite tiino Hiruims to reduce the sky signal 
in (hem before estimating the noire for a xoeund Lima, and fi- 
nally we use this improved cHtjiiuice of the uoite u> make the 
linxl maps. 

ACT muxt take a bit more care than Hll-xky experiment* be- 
C ftitw of the targe change in vcroitmiy acro^ [he map fa 
particular, near the map edgrv. only ■ « mount data 

conLitiutesi id the map, and so nrmoving the raw rap raid 
lead to an artitidal reduction m the native ol the date near 
the edge*, whicb weald lead id an artificial increase in tbor 
weight- To prevefii this from bappemitg, the muter mops used 
are first filtered and Bpodrixd. The apodizatHUi u generated 
by nrat rescaling the tm ooum map lo (be 0- 1 ,0 range, then 
vetting all values higher lhan a threshold to unity, and finally 
MikHHlung the resulting weight map with a Garovmn window. 
Furthertnore, from the Marler data map wc filter out scales 
larger Mum / ■ 300 where mmocphcric' notvc is very large, and 
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stales smaller than ( ^ 3000 where the map is highly noist- 
dom muled and simulations show tha: the bias effect is truly 
negligible f r ,t I , , we generate the starter map by multiply- 
i ng the the filtered data map with the apodiiation window. To 
test this procedure, we carry out the twostep procedure on 
both the ten] data, and the real data with a simulated map in- 
jected min it. We compare the difference of these two sets of 
maps to the original sitmilatiotu and snow the resulting, trans- 
fer function in Figure ID- At the Eow-pass scale of t = 300, the 
maps are unbiased to 0,5%, and rapidly become unbiased to 
better than a pan in I D 1 tin matter scales, 

1 1 .5. Mop Analysts W Routes 

The resulting maps cuter ami area of ^5 6 deg’ on the ’Ay, 
ranging between 20^43"* and 7 b 53 3 in right ascension and 
] and -57 2 in declination. A total of five maps were 
made 1 me using the full dataset, plus four partial versions us- 
ing independent subsets of the data for noise estimation pur- 
poses Each map is accompanied with a hit map wish the 
number of data samples thru felt on every pise I of the cor- 
responding map.. 

These maps arc an overall improvement ova: those used in 
Fowler et al. (20 10); Marriage et al, {201 La_b); Dunkley et al 
(2D1 1 ); Das et al, (2011b); Hujiun a al. (2011): Schgal el al, 
(201 U, ntariced part icn tarty by lower noise on angular scales 
between 500 < f < 2000, reducing the variance by up to a fac- 
tor of _ or 3- However, the nog mat maps were already nearly 
sample- limited on these scales, so errors on Lhe power spectre 
Lire only modestly reduced These changes are due to slightly 
di Here’ll data selection cuts, improvement in the noise treat- 
ment during the mapping process and improved pcrndetecior 
call bra uon. 



rto. IQ. — TtMivfet function of Lhe i imp (hiking pivotsc A simulucit awp 
fin) Bdiku [n the UEM itreiint, 4rtJ IlifCthw *ilh the ivM <Jsia 

(Eint-injaci imp), Then the real dun map (ptmtouity omjMud) wh ub- 
tncted from the lim-inject ttup iu pnnluu? ehewiipur mip (uutl. The master 
function wta cwppdi«t in (Ik- power spectrum of the wulpiti n«up .1 vuJi-J hy 
the crm-spKmtm of the fturpui hup end ibcirtmJiLed map. The revultlt es- 
sentially unity for all itafc* shove our nuwc-csii motion mpuL map lUtdracaiLe 
of MJ0 (>« lert), Willi « tlighl huuwtlog of large vilo that rentheu 2% It 
(nulhftfiss of 200 


Figure 1 1 shows the final 2<KJli southern map, overlaid by 
contour* of tw sensitivity (35, 5D and 65 /iK-arcmimneX as 
different region* of l he map have dt Herein integration limes. 


These sensitivities were computed using the number of sam- 
ples per unit area and the noise equivalent tempersinne given 
in Table 3. We run foe 2000 PCCi iterations on this map {o 
ensure dial it is truly converged, though as mentioned wc find 
that all science scales are converged in tar fewer PCG item- 
lions. 

Figure 12 shews the power spectrum computed for this 
maps compared to previous versions already published in 
Fowler et al (2010) and Das d al. (20! Lb), The uncertainly 
in the spectrum reveals an improve men! in the quality of the 
map 

11*2 paramaier-, obtained using u fit to the spectrum released 
with this paper are slightly shifted with respect to die values 
presented in Dunkley cl ai. (2011). The Urgesl shifts in die 
primary parameters were for n M and ft* h 1 at Acf ^0.3 and 8* 
with Ae - 0.7 t although the chApge in the Inferred dark en- 
ergy density changed by a negligible 0 03 cr, The other param- 
eters changed by less than 0.07 n . The secondary para meters 
were affected by the new mapping procedure. decreasing the 
Poisson point source level from I3,7^K 2 to 123,1 dK 1 . The 
upper 95% confidence limits ol the correlated point source 
amplitude and SZ amplitude are instead unchanged. 

Also, a cross-correlation analysis to the BLAST maps 
(htc,pi / blasi experiment . info/) show correla 
tions detected al a 25cr level implying a detection of emission 
by radio and dusty star-fotming gaieties (DSl^Gs) al high f . , 
as described in Hzjian et al. (2012). 

Finally, the dux of sources estimated in the d*ia reles*; 
is approximately 5% higher than those reported in Marriage 
rt al. (201 la), due id Lhe new source treatment in the map 
making and s change iu the beam solid angle The best rep re 
senlation of the beam shape* as well as the window function 
for power spectrum analysis, are released together with the 
maps. 


12. CONCLUSIONS 

Wc have fully character red lhe data from lhe 200ft season 
of observations from ACT. including data selection, calibre 
tkw, pot ruing, random and systematic noise and atmosphere 
contamination. 

Observations done in ZOOS yielded a total of nearly 
1260 h of CMB survey data at 148 GHz, distributed between 
two observation stripes centered at declination* near 0° and 
“53*, wilh approximately 850 and 280 square degrees respec- 
tively. After date selection, the observing time was reduced to 
nearly 8 16 h. Out of a potential 1024 live detectors, the num 
her of ell'eaivc detectors was an average 593. Combining 
observing time and detector performance, it yields a system 
efficiency of 38%. Including the ratio of observing time over 
calendar time, the overall system efficiency was 15%. 

The uncorrelated noise in live data. which excludes the at- 
mospheric fluctuations, is dominated by detector noise. Be- 
tween 5 Hz and 20 Hz, we found a total NET of 32^K v /s at 
I AS GHi. In other wools, given 1 h of observation per square 
degree, or I s of integration per square areminute, the noise 
should be about 32 }j.K areminutes- 

The curre fated noise is dominated at low frequencies by 
the almosphcre anti thermal drift, while at higher frequerp 
cfejj electromagnetic pickup ami mechanical vibrations are the 
largest sources of correlated noise. The former arc character- 
ised in frequency by negati ve power laws, equaling the detec- 
tor white noise spectra at approximately 2 Hi. The latter show 
up mostly as narrow-band signals emerging from the detector 
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Hu, I ] FiimI mop i»r i he wjuihcm reHtor observed wt MfiOB* in 34XW obtained after 2UU0 PCG iicraihms. Tt* contours in white join regions with the 
scniiK’ sensitivity, nitndy A5 h 30 nnd 65^K’iiecmlni)ti;, with inttishtitif unsiiivtiy lowuds the center at the mips The loud ruv* enclosed by the contours: is 
420 deg* end 5 MJ deg 1 respectively,. Tttt ftjiofl bHWHd bid 21*] S“ III night flsemalbn lA not shown because it was sparsely observed and hoi 

low sensitivity. Mvtl uf the anstyils tin been dura ic lIk ihqxlr itgiofl btilttiKfl roughly (1* and 7* in R,A, , Itte map tins been high-pass tillered for clarity,, 
deprttttlug USOdOA linger than t = 300- The CM II aid in j tropics Cult he sOtHt by aye. Theta is SOItie evidence Of large SOtde systematic noise, eSptriilLy near die 
edges, ii'nervad H long hurtMinlal fanturw. Which are mostly related m -hCilHiyrahn ilKnlM systematic* (HUch IK ground spillover). Despite the aijj- D f the plot, it 
it Mill possible In ih-i: Hume height point sources find SZ oIuhL-cTm 
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r-IO I J 'iJtic hy H.de coiuparistKi between <he Ad' map I '.e-uon 2-JiiH,! and ih* SPT map fur the Mtw reason of ilhe sky. The lefr [Mhef show* [he ATT nwji 
hijjh-piun littered with a ooa 3 i-ljke Eller that (oe* frttm 0 Lg I hr 1 00 C t < 300, and (lie cerdnr mid rl|{hl panels itow the ACT md SFT map* Rlpeettoiy under 
[he same higli-pa.it titter used in the 5FT d je- ;l k jic • 5>i h.r r- e< .i . 70] J ). Agreement between Ifca CMB feature! Ln the Ito { f clear by eyt 
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noijte. They can all fee modeled aa correlated modes added 10 
each elector TOD- Once identified they can be do- weighted 
in the map solution. inmimreinf! the enor> in ihe map, vtcilc 
keeptr f the solution unbiased 

Ihe dsia proeexsiog step* for mop making can be vumma- 
i laud cs ttlovt: The median ugna! kvd of each detector 11 
remov:d together with a single array -wide dope acmes the 
;5 mreute file, and the inverse anti a hi v fiber and lime con- 
stant Lzconvnhitirm me applied (sec |5 and |9). Thai ■ pre- 
calculated data selection, calibrati cm and point iris solution are 
applied (sec §10. |7 and P) After (him 1 pro computed dark 
modes are tk- projected from each detector TOD (see §1 1 and 
Appctidik A). Before map making. l be expected signal from 
the brjjht point sources U subtracted from the time streams. 
Then the maps art: made by mini mi ring \ 2 over a noise model 
which includes correlated modal, cuts (musing time streams) 
jmi frequency dependence. as described m fl M This is done 
twice, with in climate uf the map removed from the data be 
fore th; tiui*e estimation* are repealed m the second mapping 
pen. finally, the Hus, from pmnt sm.rce* previously removed 
is re -Added to the map' The retaking map* cover 845.6 dejr 
on the sky and art coflstaaN with W' MAP at angular scales, 
measured in mmnwo 

The same area of the sky covered by thew map* was ob- 
served by the South Pole Tfclesctipr (SFTi leant. who have 
recently made their rintu public l Schaffer el at, 2011}. Fig- 
ure 13 provides a aide by side comparison between the AC1' 
i i ■ 1 1 . i d ihe .SPT map for the Name region of the sky and us- 
ing the same filtering used In the SPT release, The correla- 
tion between CMB feature* in the two map* is dear to the 
eye, supporting the excellent quality of both measurements, 
ft is also clear that the ACL map is ru>h icr than the SPT map, 
which will improve in future map releases when tLna iron rhe 
folbwirtg observing rcituttu (from year* 2009 and 20 1 Oj are 
included 

Sky n-'p' for the I4>i CdEf At 1 '-■Luhero nbserv j-luos from 
2008. Jc&cnbed m this paper. are aVwUbie through NASA's 
Legacy Archive fbf Microwave Background Data Analym 


(LAMBDA), where a variety of ACT analysis software, data 
producer ami model ] cm plate* are also available. Future data 
releases will include ACT observations at higher frequencies 
and ttibscqucM ofwerviDg reasons, a- well is sky coverage in 
ACT* equatorial stripe which overlaps numerous other ob- 
serving programs 
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APPENDIX 


MODE SELECTION AND RFMOVAU 

Curt junmant -lignalsj. like the atmoaphere and sysieniaLie signaL, produce correlations belwnn dctcetorTOD* Removing Uin 
"corre'atcd" noise is important for both map making and data charactenzation, *o wc devote this appendix lo explain how this is 
done in more detail. 

A Li ncunmoM ugnat can be riHxklcd a* I lion: stream supenroptned on the TODv from a ret of individual detector* 'common 
model We call this lime stream a ■correlated mode One way to estimate it i* using an appropn^ie Iukot combmntioii of 
JWlUbJT TOO*. For instance, life common nude, defined m a linear combination with equal weight of all detector TO IX i* a 
food estimator of the atmosphere signal Organizing oil the detector TOIX mio cphireus of on * — 7 j . by nutna A. the 

currelrted mode m can be expressed as 

A#C l . <AI) 

where I ii a vector of unit magnitude representing the i^hht coeffidents of the linear combination, and s nonnali/e* the mode 
so n also has unit tringmludc Wc hencefortli ikfutc the correlated modes such that they are always normalized 

Oftet a correlated nxxie has been ide uti li rd , it can be filled and subtracted from the doia; 

A' » (l-*m r )A = A(I-« r A r Ar 1 ), (A2) 

Wc enn also construct corrclnlcd modet, ift in Other ways besides the linear combinations in Eq. A l (for example, Irom it 
ihciftK metry time stream), and in this case the iir*i cquuliiy in Ei\. A2 muai he used directly, mihcrthan the second eqiiiility. 

One ran identify important cocnhited nKdev using ihe singular value decompori Lion (SVD)of ihe matrix A: 

A = LSV\ (A3) 
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wlwre :h e cokimus of U cm be idtHtifter! us norrmlfr.ed correlated modes m, S is id diagonal matrix containing the singular values 
of A [identified as the normal faction factor, $, in Equation A1 ], and the colitm,i-Ls. of V are die eigenvectors rtf the covariance 
matrix A r A [identified as vector v in Lquaiion A I ). Hie SVD mode selection can be tailored to specific signals by first finding 
the common mode from subsets of detector TODs, and then applying the previous method to find the strongest correlated modes 
out oT iJus reduced set of modes. For cximiplc, to find modes that arc mure likely to cam-late in rows, we first find the common 
mode tram all detector TODa. in each row e' 32 row common modes in ioml|, end slack, them together as toliimm of a matrix A,** . 
fhe tf longest modes arc ilitn readily found using SVD. This method is useful to identify row and column -correlated modes, and 
sill ) truy scale mode-, troni the Himo&plicrc t or the fatter, we divide the array into 3 b square block- 1 *, finding the common mode 
in each one before applying the SVD. We call the set oF modes found this way u frmilLi-common mode" 

An import am . >ri\i deration when using modes obtained flam linear combination* of live detectors is that they are also corrc- 
lated to the sky signal. For example, naively removing the common mode effectively filters the sky, depressing scales larger than 
the army me (24* or ( < 450). Naturally, the effect is stranger for the multi common mode, which fillers scale* larger than a 
Fourth if Ehe array [6' or f . ' , tfifit)). For this reftwo. tlw multi common made is used to- calculate the drift error, tiut it cannot he 
ttfiivef] removed for making maps, nor can any oilier mode obtained from live detector*. 

In Section T I she appropriate way of projecting modes uut of the data is discussed in the context of map making, 

fittri- Mode Remm ai 

The dark detectors share the same readout as the live detectors, so they con capture systematic like thermal dnrt, cL-etromag- 
netic pickup and magnetic pickup. Moreover, ax they arc not coupled to (he optical signal, oorre fated modes obtained Iron] dark 
detector* CdLnrk modes) tan I k safely removed as u preprocessing step before making map#. For example, the ihcnnsi] drift Is well 
rcr requited by the !ow- pass -filtered common mode of the dark detectors, 

For x)w- and column *. undated electromagnctit signals, she ah; lily to identify correlated modes depends on whether ihc 
£te-.incd row nr column of ihc array contains a working djtrk detector. The army design uses one dmk detector per column, 
always in the same row. Thte ideal far identifying re I iimn -correlated modes,. but not for roweorrduted modes. For ihe fader 
we must inly on the broken live detectors with a properly working readout ciicuit* as describe,. I ul $E0, Nevenhdcts the row- 
and column correlated modes can still be identified by using SVD analysis over ail dark detectors. Nute that it is useful to first 
remove the alow thermal drift from the dank detection before trying to find these other higher frequency signals. Figure 7 shows 
an example of correlation matrices before and after having removed (he dark modes. Note dual both cohimn and row correlations 
arc significantly suppressed after the removal 
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